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Abstract: The study of compact object populations has come a long way since the determination of the 
mass of the Hulse—Taylor pulsar, and we now count on more than 150 known Galactic neutron stars and 
black hole masses, as well as another 180 objects from binary mergers detected from gravitational-waves 
by the Ligo-Virgo-KAGRA Collaboration. With a growing understanding of the variety of systems 
that host these objects, their formation, evolution and frequency, we are now in a position to evaluate 
the statistical nature of these populations, their properties, parameter correlations and long-standing 
problems, such as the maximum mass of neutron stars and the black hole lower mass gap, to a reasonable 
level of statistical significance. Here, we give an overview of the evolution and current state of the field 
and point to some of its standing issues. We focus on Galactic black holes, and offer an updated catalog 
of 35 black hole masses and orbital parameters, as well as a standardized procedure for dealing with 
uncertainties. 


Keywords: neutron stars; black holes; mass distribution; formation channels 


1. Introduction 


The recognition that extreme states of matter constitute the endpoints of Stellar Evolution 
is one of the most important achievements of the 20th century. In fact, all the concepts and 
developments of Stellar Evolution started as such in the 19th century, and evolved symbiotically 
with the new exciting “modern” Physics, General Relativity, Nuclear Physics, Quantum Me- 
chanics, and Statistical Mechanics. Stellar Evolution is their legitimate daughter and combined 
many things to create a consistent and predictive picture of how stars evolve. 

This happy confluence is particularly important for the compact remnants, leftovers of 
massive stars in which the final stages prompted matter to show its ultimate nature. This is 
how the idea of neutron stars (NS) was raised, and although the black hole (BH) concept followed 
a different path, their recognition as stellar remnants unifies the two classes as stellar corpses. 

Neutron stars are now recognized to come in many varieties. In addition to the celebrated 
pulsar group, dim isolated neutron stars (also known as XDINSs) [1-4], magnetars [5], Compact 
Central Objects in supernova remnants [6-9], Rotating Radio Transients [10,11] and some 
amazing new detections tentatively labeled as Ultra-Slow Magnetars [12] form a full family 
which we would like to understand as a whole. We shall briefly address this issue below. 

Galactic black holes, much as has been the case since the confirmation of the first one in 
Cygnus X-1, are observed generally as components of X-ray binaries, in most cases with a dwarf 
companion, although a few systems with giant companions are known. This means that most 
known BH masses have been determined from dynamical parameters (orbital period, mass 
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function, mass ratio), which have been compiled in current catalogs such as BlackCAT! [13] 
and WATCHDOG? [14]. In recent years, however, novel techniques, such as the study of quasi- 
periodic oscillations [15] and of gravitational microlensing [16,17], have allowed for masses to 
be constrained in different ways. The greatest example of this lies in the case of extragalactic 
black holes, over 100 of which have had their masses constrained from gravitational-wave 
(GW) observations of compact object mergers by the LIGO-Virgo-KAGRA Collaboration (LVK) 
3. In what follows, we will deal in detail with the current 35 well-constrained Galactic BH 
masses, and also briefly overview and compare them to the extragalactic BH mass distribution 
observed so far. 

The last twenty years or so produced in fact a large body of evidence in which the simplest 
theoretical expectations serve as an overall framework only. The process of massive star 
collapse, for example, has been deeply explored and there are now different perspectives on 
how exactly it happens, and particularly on which outcome can be expected from them [18,19]. 

At the same time, the determination of masses of NSs and BHs with good accuracy in 
larger samples has allowed a better glimpse of the astrophysical processes that lead to their 
birth. This is indeed a long-term task, since many complicated issues in the evolution of 
progenitors and explosions themselves are involved. We shall not address these issues here, but 
rather indicate some recent works that illustrate the state-of-the-art understanding of them. The 
connection with the compact star features is also dependent on the binarity of forming systems 
to a high degree, and in fact it is in binaries (in which one of the members is generally non- 
compact) where most of the masses (and some radii for NSs) have been measured. Last, but not 
least, the production of NSs by accretion induced collapse [20-22] is not properly understood 
but may be important for the whole picture, as we shall see. All these issues are entangled 
when we have to address, for example, the suggested absence of objects between 2 and 5 Mo, a 
paucity called lower mass gap in the literature [23-27]. We shall address the existence of a mass 
gap below according to analyses of the latest data [28-30]. Finally, a novel form of “seeing” 
compact objects is the now systematic monitoring of gravitational-wave events. Great insights 
on NS matter have been gathered from the event GW170817 and a substantial set of compact 
binary masses has been collected through three LVK runs, and these we also discuss briefly. 

We start in Section 2 by offering a brief overview of current issues under study with 
regard to NSs, including their demographics, maximum mass, mass distribution and radius 
measurements. We follow this in Section 3 with a more thorough discussion of the current 
state of BH observations, in particular that of Galactic objects, where we present, case-by-case, 
an updated catalog of 35 known masses, calculated in a standardized way from dynamical 
parameters. We display in Section 3.5 the resulting full Galactic BH mass distribution, and 
briefly discuss and compare it to the extragalactic distribution from LVK. Finally, we summarize 
the evolution and current state of the lower mass gap problem in Section 3.6. Our concluding 
remarks are presented in Section 4. 


2. Neutron Stars 


Neutron stars were predicted simultaneously with the discovery of the neutron itself [31], 
and related to collapsing /exploding stars by Baade and Zwicky [32] without any real proof of 
their existence. “Real” neutron stars became a reality 30 years later, when the first pulsating 
radio sources were discovered [33] and the present model (or close to it) was put forward 
by Pacini [34] and Gold [35]. However, it became clear over the years that not all NSs pulse, 
since for that rapid rotation and intense magnetic fields are necessary, according to the basic 


1 https: / /www.astro.puc.cl/BlackCAT/ 
2 https: / /sites.ualberta.ca/~btetaren/ 
3 https://www.ligo.caltech.edu/page/ligo-scientific-collaboration 
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model in which the torque is given by the electromagnetic emission and scales « B73. In 
addition, new classes of NSs have been identified, and other properties such as masses and 
radii measured with increasing precision. We shall briefly point out the main features of each 
group, the expected demographics and the inferred physical quantities in the following (for a 
recent full review, see [67]). 


2.1. Neutron Star Demographics 


It is commonplace understanding that NSs are born in massive star supernova explosions. 
However, it is almost certain that contributions from the explosions of “low-mass” massive 
stars, in the range of 8-10 Mo, which are thought to undergo electron capture onto an O-Ne-Mg 
core and leave “light” NSs is an important channel, since progenitors in this mass range are 
abundant. In addition, the rate of accretion induced collapses (AICs), either in their single- 
degenerate or double-degenerate (merger) versions, is uncertain but must be added to the total 
if some NSs are formed by them [20-22]. Simply multiplying by the Galactic lifetime the proper 
core—collapse supernova rate extrapolated from observations in galaxies similar to the Milky 
Way gives a number of NSs births of ~108 over the whole life of the Galaxy [56]. The large 
group of active pulsars today has been estimated to be around 70,000 [37,38], although this 
number quite depends on the birth parameters and evolution of magnetic fields. However, it is 
highly unlikely that the estimate is wrong by an order of magnitude, leading to the conclusion 
that most of the Galactic NSs are not pulsars, but instead belong to one of several classes that 
may be termed “hidden”, which we describe below. These are NSs that may or may not have 
functioned as pulsars in the past, but the presence of which would be much more difficult to 
establish today. Attempts to unify all NSs have been presented before [39], but as we shall see, 
new puzzling detections and unresolved problems are still ubiquitous. 

Some well-known NSs are natural candidates for the large “hidden” group. These include 
the central compact objects (CCOs) in a few supernova remnants without detected pulsations 
[6,7]. Most of the population is likely to have faded away due to its old age, well beyond the 
~ Myr scale, and therefore to establish their presence and statistics is quite difficult, although 
progress has been made; see, e.g., [8,9]. 

Detected between 1996, with the observation of RX J1856.4-3754 by [1], and 2001, when 
RBS 1774 was first observed by Zampieri et al. [2], the subpopulation termed the Magnificent 
Seven or X-ray Dim Isolated Neutron Stars (XDINSs), in addition to the isolated NS Calvera [40], 
consists of relatively close, blackbody-like cooling NSs [41] that were once considered the tip of 
the population iceberg. However, pulsations and non-zero period derivatives were detected 
over time, proving that they are actually middle-aged objects [3,4]. Nevertheless, the paucity 
of this type of NS is somewhat unexpected, and no additional candidates (except for Calvera) 
were added over many years. 

Another subpopulation which was identified unexpectedly is that of Rotating Radio 
Transients (RRATs), emitting occasionally isolated pulses of a few ms and going silent for days 
or more [11]. This narrow duty cycle makes them virtually invisible most of the time, and 
therefore their number could be very large. A simple estimate from the >100 sources already 


known is 
ins 0.5 
N =2x 10°( A ) ( ), 1 
oe 100 mJy kpc? / \ fon w 


where fon is the duty fraction and all other uncertainties have been set to unity multiplication 
factors. We see that their number could be larger than the estimated number of ordinary 
radio pulsars, but it is unlikely that the bulk of born NSs can be accommodated. In fact, the 
construction of coherent solutions for the sporadic pulses allowed an estimate of some of the 
period derivatives P, and thus an estimate of the characteristic spin-down ages. They seem 
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to be not too different from ordinary pulsars in this aspect (see Figure 1), and it has been 
conjectured that for some reason the “spark” leading to a pulse is not always operating, like 
they do in the pulsar case [10]. Of course, we do not know enough about the origin of pulses to 
discard or prove these ideas. 
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Figure 1. The P — P diagram with the new objects. The diagram had to be extended up to periods of 
~ 10° s to accommodate the ultra-long period magnetar candidates. The symbols are as follows: ordinary 
pulsars (black dots), millisecond and binary pulsars (open circles), XDINSs (green triangles), RRATs (blue 
stars), Magnetars (red squares), and the new objects in blue squares. A) GLEAM-X 162759.5—523504.3 
[12], B) AR Scorpii [47], C) J0901—4046 [46], D) J0250+5884 [42], E) J2251—3711 [43], F) J2144—3933 [44]. 


A third subgroup of importance is that of the so-called magnetars, objects in which the 
emission is related to the existence of a large magnetic field, dominating the energetics of the 
rotation (i.e., satisfying that the X-ray luminosity Ey exceeds the rotational energy IQÔ); see, 
for a recent review, [5]. Magnetars seem to possess magnetic fields a few orders of magnitude 
larger than the ordinary pulsar variety. The natural question is whether there is a continuum 
of NSs, in the sense that the magnetic fields are generated by a continuous distribution, or if 
there is some kind of gap in this quantity instead. The presence of transition pulsars, with 
magnetic fields as high as some identified magnetars seems to argue in favor of the former. 
The existence of “low-field” magnetars [45] is also suggestive of a continuum (Figure 1). In 
this way, the magnetar group was defined as the NSs in which the emission is powered by the 
magnetic field, irrespective of the precise value of its intensity. 

The latest news about the NSs population include the discovery of progressively longer 
period pulsars [46] and magnetars [12] (see caption of Figure 1). These unexpected objects are 
not only below the “classical” death line, but also challenge a definite division between the two 
groups. At some point the energy condition is employed to separate them, but otherwise their 
position in the P — P plane does not allow one to establish a clear classification. In addition, 
the first (and unique up to now) pulsar-like white dwarf AR Scorpii (object B in Figure 1) 
[47] is included in this group. For now, the latter is just a singularity, but it is likely that a 
subpopulation could be found, although unrelated to the NSs demographics. 
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What are the main questions one can ask about these subpopulations and their relation- 
ships? The list is quite extensive, and we shall only point out some of them to contribute to this 
ongoing task. 


e Do magnetic fields decay? 


This question is one of the oldest, and has been answered differently over the years. The 
idea that there must be a decay of the field within ~10° yr is quite popular, as discussed 
and modeled by Viganò et al. [48] for a diverse set of 40 sources including, e.g., the already 
mentioned magnetars, CCOs and XDINSs; see also, for a review, [49]. However, objects of this 
age still show substantial magnetic fields, and the oldest ones (“black widows”, several Gyr 
old according to evolution calculations [50]) reinforce this picture. Magnetic fields may decay 
partially to a “bottom field”, but not completely and therefore many evolutionary trajectories 
in the P — P plane proposed to explain the transit and parenthood of many subpopulations 
could be misleading. Theoretical calculations support this kind of picture [51]. 


e Which is the relationship between all these subpopulations? 


The idea that some subpopulations are just an evolutionary stage leading to some other type is 
behind all the attempts of unification [39], in one way or another. It is perhaps useful trying 
to establish which are the youngest and oldest objects. Here, we face a very general problem: 
in many cases, the age of the object is calculated with the characteristic age, and its magnetic 
field with the inversion of the electromagnetic torque equation (yielding B « VPP). However, 
it is clear by several lines of argument that this is probably an oversimplified picture. Braking 
indexes are not what the ideal electromagnetic torque predicts [52], and there is evidence for 
variations of the torque itself with time [53]. Several ideas have been put forward to link 
subpopulations; the analysis by Yoneyama et al. [54], for example, supported the idea that 
XDINSs are old magnetars, not common-type pulsars. Typical ages of magnetars (soft gamma 
repeaters [55] and anomalous X-ray pulsars [56]) are <10* yr, and some support for young 
ages comes from their associations with young supernova remnants. Other kinships have been 
suggested; for instance, Keane [57] points out that, from empirical grounds, RRATs and Fast 
Radio Bursts sources are hardly distinguishable. The clue for this association, and the possible 
relation with ordinary pulsars which are side by side with RRATs in the P — P diagram, is, of 
course, a deeper and solid understanding of the short radio pulse emission, and of ordinary 
pulsar emission in general (Kaspi and Kramer [58] present an overview of how these different 
manifestations of NSs relate to pulsars). 


e Are the new objects old magnetars? 


The ultra-slow magnetars should be, logically thinking, a latter stage in which they have cooled 
and braked to very long periods. However, their inferred magnetic fields are very high, and 
according to conventional wisdom, they should not be that old indeed. Are they actually 
related? Beniamini et al. [59] have argued that there must be a large population of ultra-slow 
magnetars in the galaxy, stressing the resiliency of magnetic fields. Is there a real difference 
between RRATs and ultra-slow magnetars? In addition, is it a complete coincidence that the 
“WD pulsar” stands nearby other confirmed NSs of this group? It is premature to give definitive 
answers to these and other important questions 


2.2. Neutron Star Mass Distribution 

The NS mass sample has 105 members now, with the addition of a few recent observations. 
An analysis of a slightly smaller sample (95 objects), presented in [29], was discussed in full by 
Rocha et al. [60], and we shall highlight the main features for completeness. 
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The first important thing is that, as agreed by several groups [61-66], the sample has 
a multimodal structure, with at least two mass scales. The analysis of [60] was performed 
assuming a Gaussian parametrization with n components, of the type 


£(mp|8) = Yr; N Compt), 2) 


where ji; and g; are the mean and standard deviation of the i-th component N, and 1; is its 
relative weight, satisfying the normalization condition );7 r; = 1. 

As discussed in Horvath et al. [67], the preferred figures for both Anderson—Darling and 
Kolmogorov—Smirnov frequentist tests are the ones appearing in Table 1. 


Table 1. p-value of two hypothesis tests for three different models. 


Model u r K-S A-D Recommendation 
Unimodal 1.48 0.35 0.025 0.032 Reject 
Bimodal 1.38, 1.84 0.15, 0.35 0.971 0.990 Do not reject 
Trimodal 1.25, 1.40, 1.89 0.09, 0.14, 0.30 0.974 0.953 Do not reject 


The p-values indicate a strong rejection of a “single mass” hypothesis (labeled as “Uni- 
modal”) and hence the confirmation of a structured mass distribution. It should be noted that 
a peak at 1.25 Mo is expected, related to the NSs produced in electron-capture supernovae, but 
the analysis does not reveal clear evidence for it. The low-mass NSs appear equally likely to 
be part of the tail of the strong maximum at 1.38 Mo, which is narrower than the second one 
around 1.8 Mo. Massive NSs belong to the latter fully. 

To visualize the form of the mass distribution, we have drawn 1000 posterior samples from 
the master sample discussed in [67], and taken their mean to compare with the maximum a 
posteriori distribution obtained. The result confirms the presence of two maxima as explained. 

A Bayesian analysis was implemented to cross-check these results. Again, the Bayesian 
likelihood is much higher for the bimodal distribution, and the mean and standard deviation of 
the two peaks quite similar, yı = 1.351 Mo, 2 = 1.756 Mo and 7 = 0.087 Mo, 02 = 0.286 Mo. 
This reinforces the conclusion that there are at least two peaks, and probably a third one 
“blended” with the objects at ~1.38 Mo, but not a single mass (called sometimes “canonical” in 
the literature, a name that is now not recommended). 

It is much more difficult to attach definite formation events to these Gaussian peaks, 
although in the long run it will be a rewarding task. The maximum mass achievable by a real 
NS, Mmax, Which has as its upper bound the Rhoades-Ruffini limit of ~3.2 Mo [68], should be 
obtained in a large sample having a statistical upper mass boundary max, since in this case, 
and provided the sample is not biased, max —> Mmax. At the present time, we can say that 
Mmax Was studied using the current sample and turns out to be ~2.5 Mo both by a simpler 37 
estimate of the second peak (Figure 2), and also by a Bayesian approach, both with and without 
an upper truncation mass introduced as an independent quantity [60]. 
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Figure 2. Gray lines represent 1000 posterior samples drawn from the sample in [29,67]. The blue curve is 
the posterior mean of these synthetic samples and the black line is the maximum a posteriori distribution. 


On the one hand, such a large value of Mmax is at odds with many estimates based on the 
analysis of the merged object in GW170817, which yield values near 2.1-2.2 Mo [69-72]. On the 
other hand, [73] have also found from GW170817 that Mmax might be as high as 2.431 0.06 Mo, 
even if the remnant is born as a uniformly-rotating NS. Rocha et al. [60] aim to avoid limitations 
arising from estimates based on single sources by obtaining their Mmax from a Bayesian analysis 
of the full sample of Galactic NS masses, which makes their resulting large mass robust. There 
are also several candidates that would be above the typical ~2.2 Mo inferred value, and one 
definite quite reliable determination of 2.35 + 0.17 Mo [74] which challenges a low Mmax 
limit and has already pushed more recent works to allow for a higher Mmax than previously 
conducted [75]. A high Mmax makes room for the lighter object in the merger event GW190814 
in the NS group, possibly at the highest achievable value of any NS mass. 

It should also be noted that, even if the statistical approach of Rocha et al. [60] is considered 
more robust than a single-object analysis, the current extragalactic NS mass sample from GW 
observations is still too small for statistics to be performed confidently upon it, with only two 
NS-NS and four BH-NS mergers [76]. Future runs of both current and future GW observatories, 
with continuously greater sensitivities, should increase the size of this sample, and eventually 
allow for an analogous procedure to be performed upon it. We point out that an alternate 
approach is to treat both the extragalactic NS and BH mass sample simultaneously with 
Bayesian methods. This is the approach of Farah et al. [28], who take into account NS-NS, NS- 
BH and BH-BH mergers, finding a break in the compact mass object distribution at 2.47)? Mo 
(3c); and the later work by Ye and Fishbach [30], who from the 4 NS-BH mergers find a lower 
limit of Mmax > 2.53 Mo with 90% confidence. The simultaneous study of both NS and BH 
masses is discussed also with regard to the lower mass gap problem in Section 3.6. 

A final remark is related to the recent announcement of a very low mass value for the 
compact object in the supernova remnant HESS J1731-347, determined using X-ray spectroscopy 


and GAIA astrometry. The value is just 0.77 +07? Mo. The reported radius, on the other hand, 


—0.17 
is small but not overly so, at 10.4 hee km [77]. The main problem is that the smallest iron cores, 


originating the lightest neutron stars, are always heavier than 1-1.1 Mo. Therefore, a much 
smaller mass would be almost impossible to accommodate unless a large fraction of the core is 
blown away in the very process of the explosion. Alternatively, this could be evidence for a 
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strange star [78], discussed and worked out for almost 40 years [79]. The difference with the 
two lowest NSs reliably determined is large, since the latter seem to have a similar mass of 
1.17 Mo [80,81]. In summary, both ends of the NS mass distribution are very interesting and 
display novelties. 


2.3. Neutron Star Radii 


After dreaming of a simultaneous determination of NSs masses and radii (see reviews by 
[63,82]), and many works in which this kind of determination was attempted indirectly, with 
uncertain results, NASA’s Neutron Star Interior Composition Explorer (NICER)* succeeded in 


measuring directly the masses and radii of two neutron stars, rendering 2.072 a ee Mo and 


12.39 = km for the massive pulsar JO740+6620 [83]; and 1.34 op Mo and 12.71 ae km for 
the millisecond pulsar J0030+0451 [84]. Even at the 1c level, the measurement of essentially 
the same radius for two neutron stars that are very different in mass means that the M — R 
sequence varies sharply around this radius, and therefore that the equation of state must be 
very stiff. The heavier the reported masses (the current record being 2.35 + 0.15 Mo reported by 
Romani et al. [74]), the stiffer equations of state need to be. In fact, the consideration of “exotic” 
equations of state complying with high masses and ~12 km radii is possible and deserves 
attention [85]. 


3. Black Holes 


The basic concept of a black hole is older than might be imagined: in the late 18th century, 
John Michell and Piere-Simon Laplace already, independently, considered the possibility of 
an object so dense that its escape velocity would exceed that of light. Not surprisingly, at the 
time, the proposal did not leave the level of pure speculation. Such a step would have to wait 
for about another 130 years, until Albert Einstein’s theory of general relativity [86] and Karl 
Schwarzschild’s solution [87] for a non-rotating, spherically symmetric, mass, which to this 
day is the basic description for the spacetime around a non-rotating and electrically uncharged 
black hole. By 1965, Ezra Newman had arrived at the general solution for a rotating, charged 
black hole, which we now call the Kerr-Newman solution [88,89]. At this time, it was already 
understood that black holes are relatively simple objects, fully defined by nothing more than 
their mass, angular momentum and charge, a statement now called the No-hair theorem. 

Although a very simple kind of object from a physical point of view, the astrophysical 
nature of a BH is one of the most complex subjects in the area. On the one hand, the fact 
that a BH can be described by only three numbers also means that we can obtain almost no 
information on the evolution of its progenitor, as it is forever lost beyond the event horizon 
[90-92]. On the other hand, the physics behind their disks, jets and magnetic fields are yet on 
the frontier of our knowledge. 

A few years after it was first predicted that collapsing giant stars could be stabilized by 
degeneracy pressure to form neutron stars (Section 2), Robert Oppenheimer and George Volkoff 
[93], in 1939, starting from the work of Richard Tolman [94], determined that even neutron 
stars have a maximum mass (the Tolman—Oppenheimer-Volkoff, or TOV, mass), beyond which 
nothing would be able to stop the collapse. Although this was a step in the direction of 
understanding black holes as one of the endpoints of stellar evolution, at the time, it was 
posited that yet another mechanism for stopping the collapse should exist. Quasars, which we 
today know to contain supermassive black holes, were first detected in the 1950s, and, in 1964, 
Yakov Zeldovich [95] and Edwin Salpeter [96] independently proposed exactly that they were 
powered by such objects, but the idea was not taken very seriously then. 


4 https: / /www.nasa.gov /nicer/ 
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In 1972, it was found independently by Thomas Bolton [97], and by Louise Webster and 
Paul Murdin [98] that the X-ray source Cygnus X-1, discovered in 1964, had a massive stellar 
companion, and from its motion a first estimate for the mass of Cygnus X-1 was derived, 
exceeding the maximum mass of a NS and making it the first stellar black hole candidate. This 
finding helped to finally convince the scientific community of the existence of black holes in the 
Universe, and since then the set of known BHs has slowly grown, along with the set of known 
BH masses, in most cases still determined from the observation of X-ray sources. Along the 
way, a rich zoo of X-ray binaries has developed (see Figure 3 for the evolution of the number of 
discovered BH X-ray transients), followed by a diversification of the methods through which 
BHs can be detected. Besides X-ray binaries, BH masses have now also been constrained in 
non-interacting binaries, and for the first time from the microlensing of background light by 
a stellar black hole. Beyond the Galaxy, the GW measurements by LVK have since 2015 built 
up a catalog of compact object mergers that as of its latest iteration, the third GW Transient 
Catalog (GWTC-3) [76] contains 90 different events with well constrained masses, 83 of which 
are confirmed BH-BH mergers, totaling 166 confirmed extragalactic BH masses. Five BH-NS 
merger or merger candidates add another 5 masses to the sample. 
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Figure 3. Cumulative histogram of discovered and dynamically confirmed BH candidates in X-ray 
transients, from the BlackCAT catalog [13], last updated in December 2021. In addition, the durations of 
X-ray missions that have discovered BH candidates in outburst are indicated. 


This diversity, however, has also meant that the treatment of mass estimates has not 
always been consistent across the field, as credibility interval conventions have changed and 
the advent of Monte Carlo methods for dealing with distributed quantities has brought about 
an ease-of-use of the entire mass probability distribution, whatever its shape, and not just of 
central values and credibility intervals. 

In this section, we build a regularized catalog of stellar BH masses and orbital parameters. 
For this, we have extensively searched the current literature and collected data both already 
present and not present in existing black hole catalogs [13,14], resulting in a total of 35 objects 
for which we have recalculated all masses, where possible, based on a standardized method and 
Monte Carlo computations. In what follows, we first discuss the treatment of uncertainty when 
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it comes to the mass estimates (Section 3.1), before briefly reviewing the nature of the observed 
systems (Section 3.2), introducing the standard mass computation procedure (Section 3.3) 
and presenting them system-by-system, along with some outstanding conflicts in their study 
(Section 3.4). In Section 3.5, we fit simple distributions to the entire mass sample, and make 
a first, simple, comparison of the resulting Galactic BH mass distribution to the extragalactic 
distribution from the LVK observations. Finally, in Section 3.6, we present a short overview of 
the current state of the lower mass gap problem in light of recent observations. 


3.1. Treatment of Uncertainties 


When probing reality, no measurement of continuous quantities is exact. With some 
degree of uncertainty always present, the best that can be hoped for is a probability distribution 
(rigorously, a probability density function, or PDF) for the actual value of each observable. 
More often than not, as a consequence of the central limit theorem, the resulting distribution is 
best described by a Gaussian, which, conveniently, is fully defined by its mean (equal to its 
mode), given as the observable’s nominal value; and its standard deviation 7, which measures 
the degree of certainty with which the observable’s actual value has been constrained. 

Even though it is a very elegant way to describe data, representing the distribution of a 
physical quantity with only two values can distort the actual measurement, especially when 
uncertainties are large and the quantity of interest has been derived from other quantities with 
their own distributions. A problem that one finds when delving into the literature of compact 
objects is the lack of clarity in the definition of each measurement. Over many different authors 
and years of publication, definitions of the credibility ranges of reported results have not 
always been explicitly given, and nominal values as well are not clearly defined to report either 
the mean or mode of the corresponding distribution, which are the same only for unimodal 
symmetric distributions. 

Dynamic mass measurements of black holes can be particularly vulnerable to this, as 
they are derived from the measurement of orbital parameters such as the period, velocity 
semi-amplitude of the companion and, of particular significance, the inclination of the orbital 
plane, on which mass estimates depend as « sin i3, and which tends to make its uncertainties 
more asymmetric. Figure 4 shows the degree to which BH dynamic masses can deviate from a 
simple Gaussian due to a broad inclination range. This gives rise to the necessity of a more 
precise, yet still easily reproducible way of describing such quantities. 
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Figure 4. 10° Monte Carlo generated masses from GRO J0422+32 (histogram), along with the Gaus- 
sian distribution N(3.97,0.95) Me from [99] (black line) and our fit of an asymmetric Gaussian 
AN (4.06, 1.08, 0.71) Mo (red line). The distributions N, AN are defined in the text. 


We therefore suggest and follow in this paper the practice of always providing any 
parameters of interest, whether directly observed or calculated, in the form of best-fit probability 
distributions of their value, so that a minimum of information is lost when using the data 
elsewhere. For simplicity, we will work with only three types of distribution: a uniform 
distribution, U(x1, x2), between x; and xz; a simple Gaussian distribution, N (u, o), with mean 
u and standard deviation 7; and what we term an asymmetric Gaussian, AN (m, 01,02), defined 
as 


72N(m, 02), ifx <m, 


x~ AN(m, 01,02) = [ (3) 


aN(m,o), if x >m, 

as a distribution over a variable x, where m is the mode of the distribution (not equal to its 
mean), while co} and o we call the superior and inferior "standard deviations". Formally, 
the parameter 0) (72) is defined st. the probability of a random variable x to be drawn from 
[m,m + 04] ([m — o, m]) is ao 
(-(co, m]). 

We highlight that, in any case, the probability distribution for the mass of an individual BH 
is a reflection solely of current limitations in our knowledge of its related observable properties 
(orbital parameters, generally). Neither the asymmetric nor the simple Gaussian are taken as 
definite best-fits to the known distributions; as stated above, a simple Gaussian often results 
naturally to good approximation as a consequence of the central limit theorem. The proposed 
asymmetric Gaussian is the simplest extension of the symmetric distribution for observables 
with asymmetrical confidence intervals. 

Equipped with these distributions, we examined all BHs with available mass measure- 
ments and carefully re-estimated their masses from known orbital parameters. We also included 
some well-constrained masses as they are reported, even if orbital parameters were not avail- 
able, as discussed next. 


~ 0.34 times the probability of x to be drawn from [m, œ) 
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3.2. Nature of Observations 


As dark objects, and except in the case of GW signals from BH-BH mergers, BHs are 
observed exclusively by their interaction with other luminous sources, which in all but one 
existing detection means a binary companion. In some of these cases, transversal and radial 
velocity (RV) measurements from an observed star can be obtained and are enough to constrain 
the presence of a binary companion in order to explain the motion. If a massive companion is 
required, but none is observed, then it is most likely to be a BH. There have been so far three 
cases in which BH masses were measured in this manner for long-period giant star-black hole 
binaries, where only the giant is directly observed [100-102]. 

For closer binaries, however, another luminous source comes into play: accretion. Either 
by filling its Roche lobe or from stellar winds, the companion loses mass, part of which is 
captured by the compact component. As this matter falls toward it, it heats up from the 
conversion of gravitational potential energy and emits chiefly in the X-rays. Although these 
X-ray sources are collectively called X-ray Binaries (XRBs), they contain a large "zoo" of binary 
classes, starting from whether the compact object is an NS or a BH. In the following, we will 
briefly review some of the terminology and important observational methods for studying BHs 
in XRBs, but much of the discussion applies also to NS-hosting XRBs. 

A first important distinction to be made among XRBs is whether accretion occurs through 
an accretion disk or not; the determining factor is if the specific angular moment J of the 
infalling matter allows for direct impact onto the compact object. If matter loses energy but no 
angular momentum, then it would orbit at the circularization radius: 


J? 

GM’ (4) 
for an accretor of mass M [103]. Typically, disk formation occurs if the accretor’s effective 
radius is smaller than Reire. For white dwarfs and NSs, the "effective radius" is equal to their 
radius if there is no significant magnetic field, but is of the order of the magnetospheric radius 
otherwise. For BHs, it is the radius of the innermost stable circular orbit (ISCO). As matter is 
accreted, both energy and angular momentum loss occur through viscosity. Once the above 
condition is fulfilled, disk formation occurs as energy is dissipated by viscosity and radiation 
faster than angular momentum is redistributed throughout the disk. The loss of energy drives 
the gas towards the lowest-energy orbit with its slowly-varying angular momentum, which is a 
circular one; thus, all the accreted gas settles into a series of concentric circular orbits, forming 
a disk. As gas keeps losing energy, the only way to fall to a lower energy orbit is by losing 
angular momentum, and in the absence of external torques, this occurs mainly by transferring 
momentum outward, causing the outer parts of the disk to spiral out [104,105]. While this skips 
over many nuances of the process, the basic picture of the accretion disk is thus formed, as an 
efficient machine for lowering highly-rotating material onto the accretor, while converting the 
orbital energy into radiation, which we can observe. 

The conditions above are nearly always fulfilled for Roche lobe overflow (RLOF), as 
the lost matter carries with it the specific angular momentum of the “parent”, and thus disk 
formation follows. Accretion from a stellar wind, on the other hand, becomes relevant for 
binaries with an O- or B-type companion in a close orbit, which is in fact the case for the 
first confirmed BH, Cygnus X-1. A simple calculation can tell us that in wind accretion the 
compact object captures only a fraction ~10~ 4-107 of the mass lost by the companion, making 
it much less efficient than RLOF [106]. It is only because wind loss rates are so large (1076- 
1077 Mo yr!) that these sources can still be observable. The winds, however, carry much less 
angular momentum, and thus do not necessarily lead to disk formation; the wind-fed pulsar 
Vela X-1, for example, has shown evidence of disk formation in the past [107], but no such 
evidence has been observed so far in the known wind-fed BHs. 


Reire = 
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Most XRBs are observed in a quiescent state, in which their behavior broadly fits within the 
picture presented above. Occasionally, however, they are observed to undergo outbursts, during 
which their luminosity is greatly increased. For BH-hosting XRBs, the clearest examples are 
called soft X-ray transients (SXTs, which can also involve a NS instead), in which quiescence 
usually lasts for ~1-50 yr with luminosities of ~10°* erg s~!, while outbursts last for ~weeks- 
months with a luminosity increase to ~10*4 erg s~! [108]. The most common picture for 
explaining these outbursts is that of a disk instability [106], based on a distinction between two 
disk states: hot and high-viscosity (outburst); or cool and low-viscosity (quiescence). The state 
of the disk is determined essentially by the degree of hydrogen ionization, thus quiescence 
requires the absence of any ionization zones within the disk. This also defines the condition 
necessary for suppressing outbursts and making the system persistent: that the disk temperature 
T always exceeds a characteristic hydrogen ionization temperature Ty ~ 6500 K. Naturally, 
disks where is always true that T < Ty would also be persistent, but very faint and lacking 
outbursts entirely; thus, if they exist, there is a strong bias against their detection [103]. We 
refer the reader to the recent review by Hameury [109] for a discussion of the disk instability 
model in the context of current problems. 

The case where the disk temperature varies between being above or below Ty brings us 
back to the SXTs. In this case, the outbursts are triggered by the formation of an ionization zone 
somewhere within the disk, which will be in the hot, viscous, state. This state of heightened 
energy dissipation favors an increased accretion rate onto the compact object, and thus an 
increased X-ray luminosity. The increased X-ray luminosity brings more and more of the 
disk into the hot state, increasing even more the rate of accretion, and so on. The heightened 
accretion tends to decrease the surface density of the gas, but the central X-ray irradiation stops 
this from effectively cooling the disk, which must remain trapped in the hot state until the 
accretion rate itself drops after a considerable accretion onto the compact object. This picture 
fits well with short-period SXTs, where the radius of the disk is limited, and it can become 
completely ionized during an outburst. Long-period SXTs, on the other hand, are able to 
comport much larger disks, a large portion of which can remain cold even during outbursts. 
The cold disk can act as a mass reservoir, making outbursts much longer than usual, to the point 
where some known persistent sources may actually be very long transients; a mass reservoir 
can also cause outbursts to occur in much quicker succession than expected. Disk warps also 
play a role, as they can affect the efficiency with which the X-ray source heats the full extent 
of the disk [110,111]. At the same time, a general “exhaustion” of the disk’s mass during long 
outbursts points towards very long quiescent periods as well, which implies a population of 
quiescent transients that have never been observed during an outburst [103]. 

The discussion above leads naturally to the main distinction to be made among XRBs: 
low-mass X-ray binaries (LMXBs) and high-mass X-ray binaries (HMXBs). LMXBs are XRBs 
with luminous components with typical masses of ~ 1 Mo or less, and so are restricted to 
RLOF accretion, which means that they virtually always contain an accretion disk (see the right 
side of Figure 5). Most, if not all, LMXBs are SXTs [103]. We also see it fit to mention here that 
the phenomenon of X-ray bursts (and their associated binaries, the X-ray bursters) are entirely 
distinct from the disk outbursts. Type I X-ray bursts are thermonuclear in nature, similarly 
to novae, and can only occur in NS-hosting XRBs. For the less well understood Type II X-ray 
bursts, disk instability models have been proposed, but are not favored (see [112] for a review 
of some current proposals). 

HMXBs contain massive companions (210 Mo) and generally undergo wind-fed accretion 
(see the left side of Figure 5), but can also accrete from RLOF, and form a disk. Even when they 
do have a disk, HMXBs rarely show disk instability outbursts because companion stars of type 
O or early B are themselves potent enough sources to keep the entire disk ionized (T > Ty). 
Outbursts can then still occur if the components are distant enough during the orbit, i.e., if the 
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period is longer than ~10 d [108], or for high-eccentricity orbits. These conditions are found in 
the particular class of Be X-ray binaries [113], in which the compact object accretes matter for 
the Be star’s own circumstellar disk; in this case, outbursts are regularly observed and probably 
caused by a burst of accretion near periastron (Figure 6). Although for some time only Be-NS 
binaries could be found, the first Be-BH binary was discovered in 2017 by Ribó et al. [114]. 


HMXB iew LMXB 


Porp = 5-1 days 


i a=25 Ro i 


Figure 5. Example of typical LMXB and HMXB. For LMXBs, accretion occurs always through Roche lobe 
overflow, while, in HMXBs, wind-fed accretion is also possible, sometimes simultaneously to Roche lobe 
overflow. Note that, while we focus on BH XRBs, the same configurations are found for NS XRBs. Figure 
obtained from Tauris and van den Heuvel [108] in Compact Stellar X-ray Sources (Copr. Cambridge 
University Press 2006), reproduced with permission of Cambridge University Press through PLSclear. 


Wind-fed HXMBs are generally persistent sources (classical HMXBs), with X-ray lu- 
minosities on the order of ~10* erg s~! [115]. However, a class of Supergiant Fast X-ray 
Transients (SFXTs) was discovered in 2005 [116,117], which shows an average X-ray luminosity 
of ~10% erg s~!, but bright X-ray flares lasting for a few days, composed of a series of bursts 
lasting for ~10° s and with Lx > 10*° erg s71 [118]. 
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Figure 6. Illustration of a Be X-ray binary. The compact object does not accrete for most of its eccentric 


orbit as the Be star does not fill its Roche lobe, and the period is too long for wind-fed accretion to take 
place. At periastron, however, the compact object crosses the circumstellar disk of the Be star and builds 
up an accretion disk, leading to a burst. Although the example indicates a neutron star, Be star-black hole 
X-ray binaries can also form. Figure obtained from Tauris and van den Heuvel [108] in Compact Stellar 
X-ray Sources (Copr. Cambridge University Press 2006), reproduced with permission of Cambridge 
University Press through PLSclear. 


We can now point toward the sources for the majority of the mass estimates presented here. 
All but 4 of the 35 systems considered are XRBs; among them, 26 are LMXBs and 5 HMXBs. 
For all HMXBs and 22 of the LMXBs, we have recalculated the black hole mass Mpy from 
known orbital parameters Pop, Kep, q, i and, when available, e. Observations of the luminous 
companion provides a powerful means to constrain the orbital period Po, and companion 
velocity semi-ampitude Kep. These constrain the mass function f(Mpy), which, together with 
the companion’s spectral type, constrains the mass ratio q. For LMXBs, outburst observations 
are also available, although, for 4 of the LMXBs, only a mass estimate was available, which we 
adopted as given. 

The inclination 7 is harder to pin down. If an absence of eclipses is confirmed, an upper 
limit can be established from q as indicated in the next section. If a disk is present, spectral 
methods for disk reflection can constrain i. When jets are present (in this case, sometimes the 
system is called a microquasar [119]), their inclination can be measured to a good degree of 
accuracy; however, it is not a safe assumption to take the jet inclination to be the same as the 
orbital inclination. In the general case of modeling the companion’s motion, detailed models 
can provide best-fit constraints on the inclination. 

Indirect observation of BHs in non-interacting binaries can still occur through the mon- 
itoring of the companion’s motion, and the masses of three BHs in our sample have been 
constrained in this manner. Observation of isolated BHs is also possible from the microlensing 
of background light, but only one such event, which we have also included, has been confirmed 
so far to have been caused by a BH. Low-frequency Quasi Periodic Oscillations (LF QPOs) in 
the X-ray spectrum provide yet another way to measure the mass of a black hole in a XRB. 
QPOs appear in the power spectrum as narrow peaks and evolve through outbursts, with LF 
QPOs appearing with <30 Hz (for a review of QPOs and their different classes, see [15]). XRB 
spectra can often be fitted by the superposition of a blackbody component and of a power 
law component [120], and Titarchuk and Fiorito [121] produced a BH mass-dependent model 


16 of 39 


which correlates the fitted photon index T of the power-law component to the observed LF 
QPO frequency v. This model has since also been used to measure the mass of XRBs in which 
LF QPOs are observed. 


3.3. Mass Computation from Orbital Parameters 


From all binaries, where available, we collect the orbital period P,,, and the velocity 
semi-amplitude Kep of the companion to the BH. From these quantities, plus the eccentricity e, 
the BH mass function is computed as 


Poke 
f (Mpa) = =a ey, (5) 


when considering a Keplerian orbit [122]. The error incurred in the BH mass from assuming a 
Keplerian orbit is still much smaller than that from observational uncertainties in Keplerian 
parameters, in particular in the orbital inclination, which might be affected by serious systematic 
errors not yet taken into account in the present work [123]. We therefore employ Equation (5) 
as the relation between binary orbital parameters throughout the work. 

In the following original sources for the LXMBs, we treat all systems as circularized, 
although a small but non-zero eccentricity might be allowed for them. For the HMXBs, 
eccentricities e < 0.01 are reported, except for Gaia BH1, which was measured to have a modest 
e = 0.45. While for Gaia BH1, we consider its eccentricity in Equation (5), for all other systems 
we treat the orbit as circular. Once again, any incurred errors from this approximation are 
expected to be inferior to those from observational uncertainties and possible systematics not 
taken into account at present. 

For each system, we also collect the binary inclination i and the mass ratio q = Mep/ Mpy. 
With these parameters, the BH mass can be determined as 


2 
MeH = Pa a A . (6) 
sin’ į 
We employ the above formula as a standardized computation of Mpgų for all systems 
with orbital parameters, even when it results in a difference to the Mpy reported by the cited 
sources. In most cases where significant differences arise, they do so only as a broader 17 
range; and any considerable differences in the central value still keep the original mass within 
1o of ours. The orbital period Po, is always an observed quantity, while Kep is taken as an 
observed quantity in all but one case, in which it is computed from P rb and f (Mpun). q is in 
some cases taken as reported by the respective source, while, in others, it is computed from the 
source’s quoted Mg and Mcp; in two cases, we re-estimate it from the companion’s spectral 
type. The inclination i is taken as reported from the sources. In all cases where the sources 
report a nominal value for a given quantity, that quantity is treated as an asymmetric Gaussian 
(Section 3.1) with mode at the nominal value; if the upper and lower uncertainties are reported, 
they are taken as 7; and 02; otherwise, these are set within 10% of the mode. When only upper 
and lower bounds are reported, the quantity is taken as being distributed uniformly between 
them. 
Whenever the inclination has not receive an upper limit while the mass ratio has been 
constrained, a lack of eclipses allows the determination of a maximum inclination from 


q 1/3 

Computations of Equations (5) and (6) are performed via Monte Carlo, and we fit our 
results to an asymmetric Gaussian, which in some cases results in a symmetric distribution 
regardless. The orbital parameters and resulting mass distributions are reported in Tables 
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2 and 3. Below, we offer a case-by-case, non-exhaustive, discussion of the objects reported 
in Table 3 and their references, under the considerations made in Section 3.1. We report the 
results as given by the references; in some cases, the mass ratio is reported as Q = Mgy/Mop, 
or the masses Mgy, Mcp themselves are given. In those cases, we use MC to convert their 
distributions to a distribution of q, unless stated otherwise. 


3.4. Collected Systems 


We present below the full catalog of BH mass measurements from 26 LMXBs, 5 HMXBs, 3 
non-interacting binaries and 1 isolated BH. We discuss each object in turn, briefly indicating 
the manner of measurement, the original results we have adopted for our catalog, and conflicts 
found during their collection. The full sample is displayed in Table 3 with sources, our 
standardized mass determination and all orbital parameters except for the eccentricity, which 
was only explicitly available for six systems; for those we display, the eccentricity in Table 2. 
We show in Figure 7 the resulting mass distributions for all 35 BHs. 


Table 2. Masses and Binary Parameters, where applicable, for six Black Hole Candidates. 


Name e References 
Cyg X-1 N(0.018, 0.003) [125] 
LMC X-1 N(0.0256, 0.0066) [124] 
M33 X-7 N(0.0185, 0.0077) [126] 
SS 433 N(0.050, 0.007) [127] 
LB-1 comp. N(0.03, 0.01) [100] 
2MASS J05215658+4359220 comp. N(0.0048, 0.0026) [101] 
Gaia BH1 N (0.451, 0.005) [102] 


Eccentricities for the seven objects for which they were given, alongside the reference. N (p, 0) is a normal distribution 
with mean ji and standard deviation g, as defined in Section 3.1. 


3.4.1. LMXBs 
Sources with orbital parameters 


We list here the 22 LMXBs for which the necessary orbital parameters for the procedure 
described in Section 3.3 were available. 


4U 1453-475 
Orosz [128] compile an inventory of reliable BH-hosting binaries at the time of writing and 
their relevant parameters. We adopted their reported P rb = 1.116407(3) d and i = 20.7 + 1.5°. 


We convert their Q = 3.2-4.0 range to a uniform q distribution. We work backwards to compute 
Kop from their f (Mgr) = 0.25 + 0.01 Mo and Pop. 


GRS 1915+105 

Greiner et al. [129] use spectroscopic data to measure Po, = 33.5 + 1.5 d and Key = 
140+15kms~!. The authors report that the jet angle determination of 70 + 2° has been 
observed to be stable over several years, so that it is reasonable to assume it as the orbital plane 
inclination i. With this inclination and a companion mass estimate of Mep = 1.2 + 0.2 Mo, they 
find Mpy = 14+ 4 Mo. We adopt their Pop, Kep and i, and compute q from their M estimates. 
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Figure 7. Violin mass plot of all 35 BHs in our sample, ordered by their mean mass (left-right, top-down). 
The horizontal lines within each of the plots indicate the quartiles of the respective distribution. 


GS 1354-64 


Casares et al. [130] present a refined treatment of GS 1354-64 (BW Cir) in relation to the first 
evidence of the presence of a BH presented in [131], and with spectroscopic and photometric 
data find Porp = 2.54451(8) d, Kep = 279.04 4.7 kms! and q = 0.121823, which we adopt. 
They report only an upper limit of i < 79° for the inclination, so we adopt U (50, 80)° from [25]. 
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GRS 1124-684 


Wu et al. [132] present a study of optical spectroscopic and photometric data of GRS 1124- 
683 (Nova Muscae 1991), and from the radial velocity of the companion are able to determine 
Porb = 0.43260249(9) d, Kep = 406.8 + 2.7 km s! and q = 0.079 + 0.007, which we adopt. We 
adopt the inclination from their following study of the object [133], in which they constrain it 
to i = (43.2777). 


XTE J1118+480 


González Hernandez et al. [134] use NIR spectroscopic data of XTE J1118+480 to mea- 
sure Pop = 0.1699337(2) d, yielding values compatible with [135,136]; and Kep = 710.0 + 
2.6 kms, also compatible with previous work [137]. We adopt their results, and their quoted 
q = 0.024 + 0.009 [138] and i = 73.5 + 5.5° [139] from previous works. 


3A 0620-003 


González Hernandez et al. [134] recalculate the orbital period derivative for 3A 0620- 
003 with new spectroscopic data from [140] and obtain a more precise value of Pon = 
0.32301415(7) d for the period, which is compatible with previous work [141]. We adopt 
their quoted Kep = 435.4 + 0.5 km s71, q = 0.060 + 0.004 [142] andi = 51.0  0.9° [143]. 


GS 2000+251 


Ioannou et al. [144] use I- and R-band photometric data to constrain the orbital parameters 
of GS 2000+251, and we adopt their nominal result of Pop = 0.344086(2) d, Q = 24+? and the 
quoted Kep = 519.5 + 5.1 km s71 from [145]. For the inclination, from their lower and upper 
limits, we adopt a uniform U(54,60)° distribution. 


MAXI J1659-152 


Torres et al. [146] use spectroscopic data of the MAXI J1659-152 X-ray transient’s quiescent 
counterpart to constrain its orbital properties and confirm the compact object’s BH nature. We 
adopt their Porb = 2.414 + 0.005 h and Kep = 750 + 80 km s71. We adopt their report lower 
and upper limits for the mass ratio and inclination as distributions U (0.02, 0.07) and U(70,80)°, 
respectively. 


MAXI J1305-704 


Mata Sanchez et al. [147] present photometric and spectroscopic data of the quiescent state 
of MAXI J1305-704 to confirm the presence of a BH in the binary. We adopt their determined 
Porb = 0.394 + 0.004 d and Kep = 5544+ 8 km s71, For the mass ratio, the authors adopt a 
N(0.05,0.02) distribution truncated to 0.01 < q < 0.07. We adapt this distribution to an 
asymmetric normal by keeping the central value at 0.05 and considering its distance to the 
truncation limits as 3g, resulting in a distribution AN (0.05, 0.007, 0.013). We keep the authors’ 
favored inclination model, i = (ay 


GS 20234338 


Casares and Charles [148] use spectroscopic data to constrain the system parameters of 
GS 2023+338. We adopt from their Table 2 Porp = 6.4714(1) d, Kep = 208.5 + 0.7 km s band 
= O0607 5 pre. Based on the parameters from [148,149], use new NIR spectroscopic data of GS 


2023+338 to measure the inclination of the system as i = (67+$)°, which we adopt. 
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XTE J1650-500 


Orosz et al. [150] use R-band photometric data from XTE J1650-500 and confirm an orbital 
period Pop = 0.3205(7) d. Reanalysis of archival spectroscopic data yields the velocity semi- 
amplitude Kep = 435 + 30 km s!. From the light curve, the authors find an inclination lower 
limit of 50 + 3° and report an upper lower limit of ~ 80° from the absence of eclipses. Although 
the exact upper limit depends on the mass ratio, we adopt a U(50,80)° distribution for i. The 
mass ratio can be measured from the companion’s spectral type, but due to the small amount 
of template spectra available, cannot be well-constrained. From the sample available, the best 
match is reported to be K4 V, from which we derive a conservative distribution of U(0.01, 0.5) 
for q. 


GRO J04224+32 


Webb et al. [151] use photometric and spectroscopic data of GRO J0422+32 to model the 
binary. We adopt their measured Kep = 378 + 16 km s7! and Pop = 0.2121600(2) d, and also 
the derived Q = 9.0*32. They constrain the mass of the compact object to >2.2 Mo, at that 
point still allowing for a massive NS. Ref. [99] studied new optical and IR photometric data of 
the X-ray transient and were able to show that it should contain a light BH, not an NS. They 
measure the most likely value for the inclination as i = 45 + 2°, which corresponds to a most 
likely Mpy = 3.97 + 0.95 Mo. While [151] find M4 V as the best-matching spectral type for the 
companion, Ref. [99] arrive at M1 V; they argue that an M4 V type cannot be made compatible 
with the observations unless an extra source of blue luminosity is posited, but that this source 
is also incompatible with observations. 


H 1705-250 


Remillard et al. [152] use photometric and spectroscopic data of H 1705-250 (Nova Ophi- 
uchi 1977) to support the presence of a BH in the system and measure P,» = 0.5228(44) d 
and Kep = 441 +6 km s71. The companion’s spectrum is best matched by a type K5, yielding 

= 0.0141) 03, which we adopt. 

Harlaftis et al. [153] use spectroscopic data from Keck to model H 1705-250 and fit also a 
K5 spectrum to the companion. They constrain the inclination to 60° < i < 80°, and we adopt 
this result as a distribution U (60, 80)°. 


GRO J1655-40 


Hernandez et al. [154] use new optical and UV, as well as archival NIR, spectrographic data 
to study GRO J1655-40 (Nova Scorpii 1994) with the main goal of performing an abundance 
analysis of its secondary. They provided updated measurements of the orbital parameters 
which we adopt: Pop = 2.62120(14) d, Kep = 226.1 + 0.8 km stand q = 0.329 + 0.047. For 
the inclination, they adopt i = 68.65 + 1.5° from [155], and so we also adopt this result. 

The authors suggest that the q = 0.419 + 0.028 found by Shahbaz [156] might be more 
accurate by virtue of having incorporated the secondary’s Roche geometry in their spectral 
analysis. However, Shahbaz [156] use i = 70.2 + 1.9° from Greene et al. [157] in their study. We 
choose to keep both the q and i from Hernandez et al. [154] for consistency. 


XTE J1859+226 

Corral-Santana et al. [158] use spectroscopic and photometric data of XTE J1859+226 and 
find from the secondary’s motion Porp = 6.58 + 0.05 h and Kep = 541 + 70 km s71, which we 
adopt. These values imply a mass function f (Mpu) = 4.5 + 0.6 Mo, in considerable excess of 
the previously reported f(Mpy) = 7.4 + 1.1 Mo [159], but we note that the older result has 
only been presented in an IAU Circular and has been considered unreliable [24]. The authors 
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also report a best-match K5 V for the companion’s spectral type, and we set for q a distribution 
U(0.01,0.5). 

The absence of eclipses imposes an upper limit i < 70° on the inclination, but the available 
data are considered not accurate enough to properly estimate the inclination angle. By assuming 
a secondary similar to that of 3A 0620-00 and the known orbital parameters, the authors arrive 
at a preferred range 60° < i < 70°. We thus adopt for the inclination a distribution U(60,70)°. 


MAXI J1803-298 


Sanchez et al. [160] use optical photometric data taken during the discovery outburst of 
MAXI J1803-298 to model the binary and provide evidence supporting the compact object’s 
BH nature. Its orbital period has so far only been measured as Porp œ~ 0.29 d, so we adopt 
a N(0.29,0.03) d distribution. The authors are able to constrain the companion’s velocity 
semi-amplitude to a Kep ~460-570 km s71 (410-620 km s7!) 68% (95%) confidence interval, so 
we take N(515,55) kms‘ as our Kep distribution. For the mass ratio, they assume a typical 
0.01 < q < 0.2 range, which we take as a U (0.01, 0.2) distribution. The authors report a >65° 
lower inclination limit but no upper limit; as no eclipses are reported, we compute an upper 
limit of ~ 85° and take U(65, 85)° as the i distribution. 


MAXI J1820+070 


Torres et al. [161] confirm the presence of a BH in the X-ray binary MAXI J1820+070 with 
spectroscopy from its decline to the quiescent state, and we adopt their Porp = 0.68549(1) d and 
Kep = 417.7 + 3.9 km s71. Ref. [162] use optical spectroscopy of the object to further constrain 
the mass ratio to q = 0.072 + 0.012 and the inclination to 66° < i < 81°, both of which we 
adopt, the inclination as U (66, 81)°. 

We note also that Atri et al. [163] provide a tighter constraint for the jet inclination as 
63 + 3°. Torres et al. [162] consider the case where the orbital inclination is taken to be the same 
as the jet inclination, obtaining a 8.48107 Mo estimate for Mgy, while their more conservative 
inclination range yields 5.73 < Mpy/Mo < 8.34. 


XTE J1550-564 


Orosz et al. [164] fit optical and NIR spectroscopic and photometric data from observations 
of XTE J1550-564 to a set of eight lightcurve models. Although the authors point out that there 
are conflicts between data taken at different times, we adopt their nominal result Pj, = 
1.5420333(24) d, Kep = 363.14 + 5.97 km s ! and i = 74.7 +3.8°; and, from their Q ~ 30, we 
adopt q = N (0.033, 0.003). Ref. [165] reports an inclination estimate of i = 37 + 4° which is 
inconsistent with the older estimate. They conclude that it may be the case that this object has 
a warped disk, or that the disk structure may be obscuring blueward line emission, resulting in 
lower inclination estimates from reflection modeling. 


GX 339-4 


Heida et al. [166] measure the donor star RV curve with data from VLT/X-shooter, chiefly 
NIR and optical, from which they obtain Pop = 1.7587(5) d and Kep = 219 +3 km s~}, which 
is significantly lower than the previous lower limit of 317 km s7! [167], leading also to a lower 
f(Mpn) than previously reported. The authors argue that the scenario in which their Kep is 
underestimated would imply variations of the absorption lines of the donor along the orbit, 
while their data cover most orbital phases and shows no such variations, which, among other 
factors, expanded upon in the cited work, support their measured Kep. 

From their Kep, Heida et al. [166] find q = 0.18 + 0.05 and a lower limit for the inclination 
of 37°. We adopt the quoted upper limit of 78° derived from the absence of eclipses and treat 
the inclination distribution as U(37,78)°. We adopt their nominal results for Porb, Kep and q. 
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GRS 1009-45 


Filippenko et al. [168] employ optical spectra from GRS 1009-45 (Nova Velorum 1993) 
during quiescence and measure Porp = 0.2852060(14) d and Kep = 475.4 + 5.9 km s~! from 
the companion’s RV curve. Approximate radial velocities of the compact primary from the 
Ha emission line allow for a mass ratio determination q = 0.137 + 0.015 and inclination ~78°. 
We adopt their results for Porb; Kcp and q, and by our standard treatment assume that the 
inclination is distributed as N(78,7.8)°. We highlight that the lack of eclipses imposes a [80° 
limit on inclination, but because the mass depends on sin? i only, which varies minimally for 
i > 80°, we keep N(78,7.8)° as our distribution for simplicity. 

Macias et al. [169] have reported refined measurements of q ~ 0.006 and Mpy = 8.5 Mo 
for this system. We, however, do not employ these results here as they are only available in the 
form of an abstract, and the cited source provides neither Kep nor f (Mgy). 


LMC X-3 


With a companion of about ~ 3.63 Mo, LMC X-3 does not fit neatly into the LMXB / 
HMXB dichotomy we have described before; it can be placed instead in the small class of 
intermediate-mass XRBs, which suffer from a negative selection bias (see [108]). Orosz et al. [170] 
study a large set of new and archival spectroscopic and photometric data to model LMC X-3. 
We employ the “Adopted value” for the parameters from the X-ray heating model in their Table 
10, Q = 1.93 + 0.20 andi = 69.24 + 0.72°. Porb = 1.7048089(11) d and Kep = 241.1 6.2 km s7! 
are taken as quoted. Their resulting Mpy = 6.98 + 0.56 Mo is considerably different from older 
estimates but is deemed more reliable as it derives from a much larger body of observations, 
and was obtained after measuring and taking into account for the first time the rotation of the 
secondary. The authors include a brief review of previous results and discuss the sources of 
inconsistencies. 


SAX J1819.3-2525 


From spectroscopic observations of SAX J1819.3-2525 (V4641 Ser), Orosz et al. [171] 
confirm the BH nature of the system’s compact object and measure Prb = 2.81730(1) d and 
Kep = 211.0+3.1 km s71, which we adopt. MacDonald et al. [172] gather 10 yr of photometric 
data from the system and are able to measure i = 72.3 + 4.1° and Q = 2.2 + 0.2, which 
we adopt. Although their estimated distance (6.2 kpc) is below that of Orosz et al. [171] 
(740 < d < 12.31kpc), their i and q were derived by fixing the same orbital period and a 
compatible Kep [173], and so we consider the parameters to be compatible. With the companion 
mass estimate as 5.49 < Mep/ Mo < 8.14 (90%) by Orosz et al. [171], this system is also an 
IMXB candidate. 


Sources with mass estimate only 


For another four LMXBs, mass estimates were available, but not their orbital parameters. 
Although this means that, for these objects, we cannot perform the standardized mass compu- 
tation from Section 3.3, we still include their masses in our sample. We treat them as symmetric 
or asymmetric Gaussians and their reported uncertainties as 17 unless stated otherwise by the 
sources. The systems here included are: GRS J1716-249, MAXI J1813-095, MAXI J1535-571, and 
XTE J1746-322. 


GRS J1716-249 

Zhang et al. [174] study X-ray spectroscopic data from the 2016-2017 outburst of binary 
GRS J1716-249 with continuum-fitting and ironline methods, with the aim of testing the Kerr 
nature of the system’s BH. For their preferred result of assuming an exact Kerr metric, they 
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find Mgy = 5.12 + 0.11 Mo, and within the same model constrain the inclination to 47.5 + 0.6°. 
We include both results in our sample. 


MAXI J1813-095 


Jana et al. [175] perform timing and spectral analyses of the 2018 outburst of binary MAXI 
J1813-095. They consider three sets of observational data and obtain a mass estimate from each 
one. We include as a mass estimate the average result, Mpy = Ale Mo (90%). They also 


place constraints on the inclination, which we adopt as a distribution U(28,45)°. 


MAXI J1535-571 


Shang et al. [176] perform timing and spectral analysis of data from its 2017-2018 outburst. 
From their spectral analysis, they find the mass range Mpy ~7.9-9.9 Mo, and we include in 
our sample their suggested mass Mpy = 8.9 + 1.0 Mo. 


XTE J1746-322 


Spectral studies of the microquasar XTE J1746-322 (H 1743-322) have been performed 
recently by Molla et al. [177] and Tursunov and Kološ [178]. Molla et al. [177] perform timing 
and spectral analyses on observations from two outbursts, in 2010 and 2011, from which 
they estimate a mass range Mpy = 9.25-12.86 Mo. They also study observed QPOs during 
the outburst, from which a second mass estimate is derived, Mga = 11.65 + 0.67 Mo. By 
combining the two results, they arrive at Mgy = 12S Mo, which we include in our 
sample. Notably, Tursunov and Kološ [178] employ a different method to also study the QPOs 
observed in this object and arrive at a mass estimate of Mpy = 11.2 Mo. 


3.4.2. HMXBs 


We list here the five HMXBs with mass estimates, all of which include the necessary orbital 
parameters for the procedure described in Section 3.3. 


Cyg X-1 

Orosz et al. [125] use an improved distance measurement for Cyg X-1 together with previ- 
ously published optical data to model the system. The period is fixed to Porp = 5.599836(24) d 
from [179] and the data are fitted to four models. We adopt the fixed period and the final 
parameters from Table 2: i = 27.06 + 0.76°, e = 0.018 + 0.003, Mgy = 14.81 + 0.98 Mo and 
Mep = 19.19 + 1.90 Mo. We compute q with MC from Mgy and Mcp. We adopt their more 
scattered Kep = 75.57 + 0.70 km s7! fitted to the data from [180]. 


LMC X-1 


Orosz et al. [124] use optical spectroscopic, and optical and NIR photometric, data to 
model LMC X-1 and confirm earlier work by Hutchings et al. [181] with much higher precision. 
We employ the “Adopted Value” from their Table 3 for Pap = 3.90917(5) d, Kep = 71.61 + 
1.10kms~!, Mpy = 10.91 + 1.54 Mo, Mep = 31.79 + 3.67 Mo and i = 36.38 + 2.02°. We also 
report their resulting e = 0.0256 + 0.0066 from their "Eccentric Orbit" values in Table 3. 


M33 X-7 


Orosz et al. [126] model the eclipsing binary M33 X-7 with optical spectroscopic and 
photometric data, fixing the orbital period to the previously determined Po, = 3.453014(20) d 
[182]. We adopt their selected parameters from Table 2: Kep = 108.9 + 5.7 km s71, e = 0.0185 + 
0.0077, i = 74.6 + 1.0°, Mgy = 15.65 + 1.45 Mo and Mep = 70.0 + 6.9 Mo. 
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NGC 300 X-1 


Crowther et al. [183] confirmed the nature of system NGC 300 X-1 as a Wolf—Rayet/black 
hole binary located in the Sculptor group galaxy NGC 300 (d = 1.88 Mpc). They use optical 
spectroscopic data and determine Porp = 32.3 + 0.2 hand Kep = 267.5 + 7.7 km s71, of which 
we adopt the latter. For the period, we adopt the more recent Pj, = 32.7921(3) h from 
Binder et al. [184], where new X-ray and UV observations are combined with archival X- 
ray observations to further constrain the binary model. We adopt their mass estimate of 
Mep = 26*7 Mo for the WR companion and Mgy = 17 + 4 Mo to recover the q distribution; 
these estimates are consistent with 60° <i < 75°, so we adopt a U(60,75) distribution for the 
inclination. 


SS 433 


Discovered in the 1970s, SS 433 is the first observed Galactic microquasar, today generally 
agreed on to be an eclipsing X-ray binary undergoing supercritical accretion onto the compact 
object. Over the more than 40 years during which this object has been studied, considerably 
different estimates of the component masses have been put forward. We have taken our 
parameters for this system from the recent review by Cherepashchuk [127], which reports 
results from the over 40 years of observations of SS433 in the optical, radio, and X-rays. 
We adopt the quoted nominal value for the orbital period from the review and assign it a 
conservative uncertainty in the last digit, yielding Po: = 13.082(5) d. The inclination is given 
as i = 79°, so we adopt the distribution N(79,7.9)°. We also collected the measured eccentricity 
of e = 0.050 + 0.007. 

For the velocity semi-amplitude, we adopt the reported Kep = 58.2 + 3.1 km s~! from 
Picchi et al. [185], and, for consistency, also adopt the mass estimates from that same work, 
Mpy = 4.2 + 0.4 Mo and Me = 11.3 + 0.6 Mo, which result in a mass ratio q ~ 2.64. 

A more recent estimate by Cherepashchuk et al. [186], however, has led to a lower limit 
of q > 0.8 which implies a relatively massive BH with Mgy > 8 Mo. We do not adopt this 
result presently as it does not provide the parameters necessary for our standardized mass 
computation. 


3.4.3. Other Sources 


In addition to the X-ray binaries, we have collected four recent BH observations with mass 
and, where applicable, orbital parameter measurements, which we detail here. These are the 
three non-interacting binaries, LB-1, 2MASS J05215658+4359220 (J0521+435) and Gaia BH1; 
and the microlensing event OGLE-2011-BLG-0462. 


LB-1 comp. 

Liu et al. [100] report the results from over two years of radial-velocity measurements 
of the Galactic B-type star LB-1 (LS V +22 25) and find that its motion requires the presence 
of a Mg = 68713 Mo (90%) BH ina Porp = 78.9(3) d orbit. From the best-fit model for the 
Ha emission line, they also find Kep = 52.8 + 0.7 km s7! and e = 0.03 + 0.01, and, from the 
companion’s spectra, they find Mcp = 8.29 Mo. We adopt these values for the object. 


2MASS J0521+435 comp. 

Thompson et al. [101] study spectroscopic and photometric data from the giant star 
J0521+435 and determine that it must be in a Po, = 83.20(6) d, Kep = 44.6 + 0.1 km s71 
and e = 0.0048 + 0.0026 orbit with and an unseen BH companion. By comparing the giant’s 
properties to single-star evolutionary models, they find for the inclination sini = (0.97+8-93)° 
and for the masses Mep = 3.2119 Mo Mpy = 3.353 Mo (20). We adopt these nominal results 
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for the object, although the authors also employ other methods that arrive at different but 
compatible estimates; the unseen companion’s mass is consistently estimated to be in the 
~2.9-4.0 Mo range, which keeps it a likely low-mass black hole. 


Gaia BH1 


El-Badry et al. [102] study a a bright solar-type star in Ophiuchus with astrometric data 
from the Gaia mission and follow-up spectroscopy, and determine from the Gaia orbital 
solution and RV measurements that the star must orbiting an unseen companion. We adopted 
the quoted Mcp = 0.93 + 0.05 Mo mass for the solar-type companion, and their results derived 
from astrometric and RV data: Por» = 185.59(5) d, Mgy = 9.62 + 0.18 Mo, i = 126.6 + 0.4° and 
e = 0.451 + 0.005. We adopt the Kep = 66.7 + 0.6 km s71 derived from RV data only. This is the 
only confirmed BH binary with a modest but non-negligible eccentricity, as well as the closest 
known BH, at the moment of writing, and potentially originated from either a triple system or 
dynamical assembly. 


Table 3. Masses and Binary Parameters, where applicable, for 35 Black Hole Systems. 
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Name Porb Kep f (Mpu) i q Megu Referenc 
day km s~! Mo deg Mo 

4U 1543-475 N (1.116407, 0.000003) AN(129.33, 1.68, 1.76) N(0.25,0.01) N(20.7,1.5) u(0.25, 0.31) AN(8.78, 2.45, 1.45 [i 
GRS 1915+105 N(33.5,1.5) N(140, 15) AN(8.79, 3.69, 2.35 N(70,2) PN 08, 0.04, 0.02) AN(12.92, 5.39, 3.57) m: 
GS 1354-64 N (2.54451, 0.00008) N (279.0, 4.7 AN (5.72, 0.30, 0.28 u (50, 80) N (0.12, 0.03, 0.04) AN (8.23, 3.76, 0.84 [25,1 
GRS 1124-684 N(0.43260249, 0.00000009) N(406.8, 2.7 N(3.018, 0.06) AN (43.2, 2.1, 2.7) ANOA. 0.007) AN(10.70, 1.98, 1.06) [132,1: 
XTE J1118+480 N(0.1699337, 0.0000002) N(710.0, 2.6 AN (6.304, 0.069, 0.07) N(73.5,5.5) N(0.024, 0.009) AN(7.23, 0.88, 0.37 [1 
3A 0620-003 N(0.32301415, 0.00000007) N(435.4,0.5 AN(2.762, 0.009, 0.01) N(51.0,0.9) N(0.060, 0.004) AN(6.60, 0.27, 0.24 [1 
GS 2000+251 N (0.344086, 0.000002) N(519.5,5.1 AN(5.00, 0.15, 0.14 u(54, 60) AN (0.034, 0.015, 0.006) AN(8.97, 0.84, 0.51 [1 
MAXI J1659-152 N (0.1006, 0.0002) N (750,80) AN(3.58,3.15, 1.91 uU (70,80) u (0.02, 0.7) AN (4.25, 3.88, 2.28 [1 
MAXI J1305-704 N (0.394, 0.004) N(554,8) AN(6.92, 0.32, 0.29 AN(72,5,8) AN (0.050, 0.007, 0.013) AN (8.18, 1.87, 0.57 [1 
GS 2023+338 N (6.4714, 0.0001) N (208.5, 0.7 N (6.078, 0.061) AN(67,3,1) AN (0.060, 0.004, 0.005) AN (8.58, 0.31, 0.45 [148,1 
XTE J1650-500 N (0.3205, 0.0007) N (435,30) AN (2.64, 0.64, 0.48 u (50, 80) u (0.01, 0.5) AN(5.03, 3.04, 1.36 [1 
GRO J0422+32 N (0.2121600, 0.0000002) N(378,16) AN(1.19,0.16, 0.14 N(45,2) AN (0.116, 0.079, 0.071) AN (4.06, 1.08, 0.71 [99,1 
H 1705-250 N(0.5228, 0.0044) N(441,6) AN(4.64, 0.20, 0.19 u (60, 80) AN (0.014, 0.019, 0.012) AN(5.41, 1.15, 0.36 [152,1 
GRO J1655-40 N (2.62120, 0.00014) N(226.1,0.8 N(3.139, 0.033) N(68.65, 1.5) N(0.329, 0.047) AN(6.86, 0.52, 0.49 [154,1 
XTE J1859+226 N(0.274, 0.002) N(541,70) AN(4.08, 2.12, 1.29 u (60,70) u(0.01,0.5) AN(7.79, 5.70, 2.64 [123,1 
MAXI J1803-298 N(0.29,0.03) N(515,55) AN(3.77, 1.66, 1.04 Uu(65, 85) u(0.01, 0.2) AN(5.06, 2.51, 1.45 [1 
MAXI J1820+070 N (0.68549, 0.00001) N(417.7,3.9 AN(5.18,0.15, 0.14 uU (66,81) N (0.072, 0.012) AN (6.47, 0.80, 0.30 [161,10 
XTE J1550-564 N(1.5420333, 0.0000024) N (363.14, 5.97) AN(7.64, 0.39, 0.36 N(74.7,3.8) N(0.033, 0.003) AN(9.00, 0.78, 0.56 [ie 
GX 339-4 N(1.7587, 0.0005) N(219,3) AN(1.913, 0.08, 0.007) u(37,78) N(0.18, 0.05) AN(3.21, 3.50, 0.44 [i 
GRS 1009-45 N(0.2852060, 0.0000014) N(475.4,5.9 N(3.17,0.12) N(78.0, 7.8) N(0.137, 0.015) AN (4.28, 0.57, 0.24 [i 
LMC X-3 N(1.7048089, 0.0000011) N(241.1,6.2 AN(2.47, 0.20, 0.18 N (69.24, 0.72) AN (0.512, 0.099, 0.084) AN(7.20, 1.23, 0.94 [uv 
SAX J1819.3-2525 N(2.81730, 0.00001) N(211.0,3.1 N(2.74, 0.12) N(72.3,4.1) AN (0.447, 0.048, 0.034) AN (6.64, 0.79, 0.54 [171,1; 
GRS 1716-249 : > N(47.5,0.37) Š N(5.12,0.11) p; 
MAXI J1813-095 - - - U(28, 45) - AN(7.41, 0.90, 0.93 [i 
MAXIJ1535-571 2 r : 4 : N(8.9,1) [1 
XTE J1746-322 - - - - - AN(11.21, 1.65, 1.96) [uv 
Cyg X-1 N(5.599836, 0.000024) N(75.57,0.70) N(0.250, 0.007) N(27.06, 0.76) N(1.29, 0.17, 0.14) AN(13.72, 2.55, 1.88) [125,17 
LMC X-1 N(3.90917, 0.00005) N(71.61, 1.10) AN(0.143, 0.005, 0.004) N (36.38, 2.02) N(2.81, 0.62, 0.41) AN(9.27, 3.85, 2.03 [1: 
M33 X-7 N(3.453014, 0.000020) N(108.9,5.7) AN (0.449, 0.082, 0.062) N(74.6,1.0) N(4.39, 0.68, 0.53) AN (14.39, 5.12, 3.13) [126,1% 
NGC 300 X-1 N(1.3663375, 0.0000125) N(267.5,7.7) AN (2.68, 0.26, 0.23) U(60,75) N(1.42,0.70,0.33) AN(19.79, 14.01, 5.38) [183,15 
SS 433 N(13.082, 0.005) N(58.2,3.1) AN (0.263, 0.047, 0.038) N(79.0, 7.9) N(2.64, 0.34, 0.24) AN (4.23, 1.70, 0.98 [127,18 
LB-1 cp. N(78.9,0.3) N(52.8, 0.7) AN(1.203, 0.049, 0.047) u(15, 18) AN(0.117, 0.022, 0.019) AN (56.86, 18.59, 4.75) [i 
2MASS J0521+435 cp. N(83.2, 0.06) N(44.6,0.1) N(0.765, 0.005) AN(75.43, 5.16, 2.35) AN(0.74, 0.26, 0.17) AN (2.45, 0.84, 0.45 [i 
Gaia BH1 N(185.59, 0.05) N(66.7,0.6) N(4.08, 0.08) N(126.6, 0.4) N(0.097, 0.005) AN(9.42, 0.33, 0.31 [i 

é 


OGLE-2011-BLG-0462 


N (7.88, 0.82) 


[ 


In each line, we identify the name of the source containing the black hole, as well the BH mass and available orbital parameters; when these are available, the indicated BH mass is the result of our standardized calcula 
(Section 3.3); otherwise, it corresponds to the nominal result from the references. For the systems with more than one source, which parameters where taken from which source are specified in Section 3.4. The distribut: 
are given in terms of the notation established in Section 3.1: N (u, 0) indicates a normal distribution with mean yi and standard deviation 7, while AN (m, 01,02) an asymmetrical normal distribution with mode m, stanc 
deviation dı above m and standard deviation oz below. U(x1, x2) indicates a uniform distribution between x1 and x2. 
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OGLE-2011-BLG-0462 


The microlensing event OGLE-2011-BLG-046S (OB110462) marks the first discovery an 
isolated BH. It was made independently by Lam et al. [16] with data from the Optical Gravita- 
tional Lensing Experiment (OGLE), and by Sahu et al. [17] with data from the Microlensing 
Observations in Astrophysics (MOA) survey, with initially conflicting mass estimates of 1.6- 
4.4 Mo and 7.1 + 1.3 Mo, respectively. A simultaneous re-analysis of all available data from 
OB110462 by Mróz et al. [187] later confirmed a mass of Mpy = 7.88 + 0.8 Mo for the object, 
which we include in our sample. 


3.5. The Full Mass Distribution 


One of the great interests in building an updated and standardized catalog of BH masses 
is to verify the shape of the full Galactic BH mass distribution and to match it to analytical 
models, which may be made available for use in further work concerning BH populations. Ozel 
et al. [24] and Farr et al. [25] performed the last such modeling more than 10 years ago, with 
samples of 16 and 20 BHs, respectively. At that time, through Bayesian methods, they found 
that, for a simple Gaussian, the preferred distribution had a mean u ~7.8 Mo and standard 
deviation 7 ~1.2 Mo. On the other hand, Farr et al. [25] also considered a double Gaussian 
model, and found for it a Bayesian evidence ~ 30 greater than that for the simple Gaussian. We 
present here a simple goodness-of-fit test for an asymmetric Gaussian model to our full sample 
and compare it to the previous results. A full investigation of the mass distribution with more 
robust methods and including further corrections for systematics and biases [123,188] will be 
developed in a future work (Bernardo et al., in preparation). 

From each of the 35 mass distributions in Table 3, we draw 10° masses, and fit an asym- 
metric Gaussian to the mixture of these draws, resulting in m = 6.72 Mo, 0 = 4.17 Mo and 
02 = 2.27 Mo. The sample histogram is displayed along with the asymmetric Gaussian fit in 
Figure 8. The mode of the distribution (m) results in a value not too distinct to the result from 
[24] for a sample two times smaller, indicating the robustness of a peak in this region, while 
allowing for asymmetry, leads to quite different c. 
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Figure 8. Mixed sample resulting from a draw of 10° masses from each object in Table 3 (yellow histogram) 
and the resulting best asymmetric Gaussian fit as AN (6.72, 4.17, 2.27) Me (blue line). 


While a full analysis of the currently known extragalactic BH masses measured by LVK is 
out of the scope of this work, some first considerations and a simple comparison stand to be 
made. It must first be highlighted that there is no simple way to derive a single distribution 
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from both Galactic and extragalactic samples simultaneously, as not only are they affected 
by different observational biases, but also the underlying physical distribution itself is not 
necessarily the same. Different environmental conditions at the higher redshifts probed by GW 
observations could change the shape of the distribution by affecting the process of stellar and 
binary evolution. Lower metallicites, for example, make wind mass loss from giant stars less 
efficient (e.g., [189]), while both low metallicities and high star formation rates have been linked 
to a top-heavy initial mass function (see the reviews by [190,191]), which could be expected to 
lead to a BH mass distribution shifted to higher masses. 

With that in mind, we have collected the raw BH mass distributions made available 
by LVK (Monte Carlo samples of 8,000 masses for each object) as part of GWTC-3 [76] and 
produced their mixture distribution with no further processing of the data. Unlike the Galactic 
distribution, the extragalactic distribution shows two peaks, which we fitted to a double 
Gaussian in order to determine the location and breadth of each one. Figure 9 shows the 
resulting fit, with a lower peak p1 = 9.07 Mo, 01 = 2.39 Mo, with weight r; = 0.17; and 
an upper peak u2 = 33.86 Mo, 02 = 16.85 Mo. Even though a direct comparison cannot be 
directly made, as discussed, it is still noteworthy that a lower peak from the GW sources seems 
to align with the single Galactic peak as shown in Figure 10. This lower peak is mainly caused 
by the new additions to GWTC-3 in relation to GWTC-2 [192]. It may then be connected to the 
improvements in sensitivity made to the observatories between each new run, which allow for 
the detection of lower-mass sources than in previous runs. Although a safe conclusion may 
only be derived after a careful modeling and consideration of observational biases for both 
samples, this might lend credence to the idea that a Galactic-like population of BHs exists in 
galaxies at higher redshift (z~1). 


Combined GW Black Holes 
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Figure 9. Mixed sample resulting from a draw of 10° masses from each BH component from GWTC-3 
(yellow histogram) and the resulting best double Gaussian fit (blue line), with parameters given in the 
text. 
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Figure 10. Comparison between the same mass draw from Galactic BHs shown in Figure 8 (blue 
histogram) and a mass draw from the extragalactic BHs from GWTC-3 (yellow histogram). We note the 
alignment between the Galactic peak and the lower extragalactic peak. 
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Figure 11. BH-BH merger rate (RBBH) per total binary mass (Mtot) bin in three different ranges of redshift 
at merger Zmerger from [193], resulting from binary population synthesis with an IMF that becomes 
top-heavy at high z. 


On the other hand, the GW mass distribution shows a high-mass peak which is entirely 
absent from the Galactic distribution. Although it is hard to treat the existing Galactic BH 
sample as reflective of the real mass distribution with its current size, this also might indicate 
that high-redshift populations allow for the production of more massive BHs than local popu- 
lations. If mass-redshift correlations are present in the GW sample, then they might support 
the aforementioned idea that low-metallicity and high-star formation rate conditions at high 
redshift lead to the birth of more massive stars, [194] which in turn allow for an increased 
production of more massive BHs, and consequently of more massive BH mergers, as has been 
indicated by population synthesis studies that consider environment-dependent initial mass 
functions and star formation rates ([193,195], see also Figure 11) and tentatively supported 


30 of 39 


by GW observations [196]. A full consideration will, however, still require accounting for 
observational biases and for a study of the redshift distribution in GWTC-3 and future catalogs. 


3.6. The Lower Mass Gap 


We note from Figures 8 and 9 that both the Galactic and extragalactic mass distributions 
smoothly extend into the <5 Mo region to a significant level; in fact, six of our objects — MAXI 
J1659—152, GRO J0422+32, GX 339—4, GRS 1009—45, SS 433 and 2MASS J0521+435 cp. — all 
have mass distributions with central values below 5 Mo. This makes them candidates for being 
within the 2-5 Mo lower mass gap in the BH mass distribution, an apparently empty mass range 
above the maximum NS mass first indicated by [23] from a sample of seven LMXBs. With their 
expanded samples of 20 XRBs, Ozel et al. [24] in 2010 and Farr et al. [25] in 2011 still found 
evidence for a gap in the BH mass distribution below ~5 Mo, and so its existence has remained 
an open question for over 20 years. The gap has been connected to the characteristic convection 
timescale in core-collapse supernovae [26,27], but even in this case a partial filling has not been 
ruled out. Ongoing work on this issue [197,198] can now be taken both ways, as while in one 
hand an appropriate convection timescale can help explain the existence of the gap, a measure 
of the gap depth can help constrain the convection timescale. Other mechanisms for compact 
object production, such as AIC in WD binaries and accretion in "spider" binaries, have also 
been suggested to fill the gap to some extent [67], but the hypothesis of a partially filled gap 
has remained mostly speculative until recently. 

Objects such as the six we mention from our sample, as well as light BHs detected by 
LVK have in the past few years started to offer observational support for the partial gap 
hypothesis. Farah et al. [28] performed a Bayesian study of the full sample of BH and NS 
masses measured by LVK and considered a single, overall compact object distribution in the 
form of a broken power law with a parametric dip representing the gap. Although they cannot 
conclusively confirm the existence or absence of the gap, their results show a preference for 
a partially depleted gap with a preferred lower edge at 2.20. Mo, not incompatible with 
~2.5 Mo maximum NS mass estimates. The broken power law showed a preference for a break 
at 2.440 Mo, which could point towards a high NS Mmax- 

de Sá et al. [29] later performed a frequentist analysis of 12 current gap candidates, both 
Galactic and extragalactic, in relation to the then most recent BH and NS mass distributions 
— which show a gap between ~2.2-5 Mo — displayed in Figure 12. They calculate both the 
likelihood of individual objects falling in the gap, and by a likelihood ratio test between the 
existing distributions and a simple partially depleted gap model are able to show that current 
detections place the probability of a desert gap below ~3 Mo as being at most ~10%, and that 
a mass distribution with a partially depleted gap is strongly preferred over one with a desert 
gap, with a likelihood ratio of ~10%. It seems thus clear at this point that, if there is a gap, it 
should not start below ~2.5 Mo, and that the investigation of alternative formation channels 
for compact objects is needed to further constrain the degree to which the 2.5-5 Mo range is 
populated. 
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Figure 12. Mass distribution and gap object candidates from [29]. The solid line shows the mixture of the 
BH mass distribution from [24] and the NS mass distribution from [60]. In addition, shown are the mass 
measurements, with their 17 confidence intervals, of the 12 gap candidates the authors consider, a mix of 
Galactic and extragalactic sources. 


4. Conclusions 


The size and density of the available compact object sample has greatly grown since 
their existence was first confirmed in the 1960s, with the sample of NSs surpassing 100 and 
that of Galactic BHs reaching 35 objects, while the GW observations from LVK have already 
detected more than 150 BHs and around 10 NSs at high redshift. Slowly, this growing sample 
has allowed us to start to see the shape of the underlying compact object mass distribution, 
provided biases do not dominate the samples. Once thought to be clustered around a canonical 
mass of ~1.4 Mo, today NSs have been shown beyond doubt to have masses >2.2 Mo, with at 
least two peaks in the mass distribution. More recently, growing evidence suggests that their 
maximum mass might be closer to ~2.5 Mo. This slow upward shift of the maximum NS mass 
value serves also as a cautionary tale with regard to the BH mass distribution, which is still 
based on only a few tens of measurements. As the canonical model for NS masses changed. 
considerably while the sample was below ~100 objects, the BH sample may still need an 
accumulation of reliable data points before it can properly constrain the shape of the physical 
mass distribution. It is expected that >10” BHs and a higher number of NSs should exist in 
the Galaxy, originating from a great diversity of evolutionary paths which we are only starting 
to understand. Progress in this front will also be needed before we can properly evaluate the 
degree of completeness of the sample due to observational biases. 
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